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I . SUMMARY 


This  report  summarizes  the  technical  results  obtained 
under  RADC  Contract  No.  F 30602-72-C-0386  during  the  oeriod 
March  1,  1974  through  March  31,  1975.  It  is  the  final 
report  on  the  University  of  Colorado's  effort  in  support 
of  the  ARPA  HF  Adaptive  Array  Program.  The  work  is  a con- 
tinuation of  that  reported  earlier  and  has  been  carried 
out  in  close  coopercition  with  the  Remote  Measurements 
Laboratory  at  Stanford  Research  Institute,  Menlo  Park, 
California.  Data  for  the  study  were  obtained  using  the 
WARF  which  is  located  in  the  central  valley  of  California 
and  which  is  operated  by  SRI. 

Results  have  been  obtained  relating  to  the  applicability 
of  time-domain  adaptation  procedures  in  HF  backscatter 
radar  systems.  Specifically,  investigations  regarding  the 
performance  of  these  procedures  in  an  actual  HF  interfer- 
ence environment  were  conducted  in  four  distinct  areas: 


1.  Coefficient  Time-Variation  Studies: 

It  was  shov/n  that  the  time  scale  of  coefficient  varia- 
tion in  an  HF  environment  was  the  order  of  about  one  second. 
These  results  indicate  that  block  adaptation  methods  such 
as  matrix  inversion  techniques  must  use  an  update  interval 
of  less  than  one-half  second.  It  was  further  shown  that 
the  P-vector  adaptation  method  can  successfully  track  such 
variations  in  a 40  weight  processor  and  that  the  standard 




i 


deviation  of  the  resulting  weights  is  accurately  predicted 
by  adaptive  theory. 

2.  Clutter  Suppression  Studies 

The  use  of  individual  three-pulse  MTI  clutter  suppression 
filters  at  each  subarray  output  has  been  investigated. 

For  the  experimental  data  studied,  no  significant  advantage 
or  disadvantage  was  observed  with  the  use  of  these  filters. 
Typical  observed  variations  were  the  order  of  a few  dB. 

3.  Aperture  Tradeoff  Studies 

It  was  shown,  using  a simple  simulation  model  and  exper- 
imental data,  that  in  a non-isotropic  noise  environment, 
conventional  processing  does  not  provide  uniformly  in- 
creasing performance  with  increased  aperture  size.  It  was 
further  demonstrated  that  with  adaptive  processing,  the 
improvement  is  a monotonically  increasing  function  of 
aperture.  The  presence  and  effect  of  fading  nulls  on 
aperture  utility  was  observed  using  experimental  data. 

Such  nulls  can  significantly  affect  the  selection  of  optimal 
subarray  location  and  spacing.  Several  examples  demonstrat- 
ing this  phenomenon  are  presented. 

4.  Minimum  Mean-Square  Error  Adaptation  with  Constraints 
An  alternative  time-domain  adaptation  method,  suggested 

by  Frost,  has  been  implemented  md  compared  with  the  P-vector 
algorithm.  Frost's  procedure  has  the  specific  advantage 
of  providing  a flat  frequency  response  in  the  main  lobe 
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direction. 


It  is  demonstrated,  however,  that  the  extra 
degrees  of  freedom  required  to  produce  this  response 
result  in  a lower  output  signal-to-total-noise  ratio  than 
that  provided  by  P-vector  adaptation.  The  amount  of  this 
degradation  was  observed  to  be  as  much  as  ten  dB  in  some 
range-doppler  cells. 

Complete  details  of  these  results  are  presented  in  the 
sections  following. 
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II.  INTRODUCTION 


This  report  is  the  final  report  on  the  University  of 
Colorado's  effort  in  support  of  the  ARPA  HF  Adaptive  Array 
Program.  The  results  described  herein  represent  a continuation 
of  that  described  in  References  1,  2,  and  3.  These  references 
contain  complete  descriptions  of  the  experimental  equipment 
used  for  data  collection  as  well  as  mathematical  details  re- 
garding the  structure  of  the  P-vector  adaptation  algorithm 
used  throughout  this  study. 

The  overall  objective  of  the  HF  Adaptive  Array  Program  has 
been  to  investigate  the  degree  to  which  adaptive  beamforming 
methods  can  improve  the  signal-to-noise  ratio  (SNR)  observed 
at  the  output  of  an  HF  backscatter  radar  system.  Efforts  at 
the  University  of  Colorado  have  been  directed  toward  the  devel- 
opment and  analysis  of  specific  adaptive  algorithmic  techniques 
designed  to  meet  this  objective.  Data  for  the  study  were  re- 
corded using  the  Wide  Aperture  Research  Facility  (WARF)  which 
is  operated  by  the  Stanford  Research  Institute,  Menlo  Park, 
California.  Personnel  at  SRI  have  worked  in  close  cooperation 
with  the  University  of  Colorado  on  this  program  and  have  devel- 


oped and  tested  an  on-line  adaptive  processing  system  for  use 

[4] 


at  WARF.  A recent  SRI  report  L"'J  summarizes  their  effort. 

The  results  summarized  in  the  present  report  represent 
further  investigations  into  the  use  of  time-domain  adaptation 
procedures  in  HF  array  applications  and  include  the  following 
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specific  studies: 

1.  A study  of  the  nature  of  the  time-dependent  behavior  of 
adaptive  weighting  coefficients  in  the  processor  during 
adaptation  on  field  recorded  data. 

2.  An  investigation  into  the  use  of  clutter  suppression  filters 
(MTI  filters)  at  the  subarray  receiver  outputs  and  prior 

to  adaptation. 

3.  An  aperture  study  designed  to  measure  performance  versus 
receiving  array  aperture  with  an  adaptive  beamformer  and 
to  compare  these  results  with  those  obtained  using  con- 
ventional beamforming  methods. 

4.  A comparison  of  the  P-vector  adaptation  algorithm  with  a 
constrained  least-squares  adaptive  procedure  suggested  by 
Frost.  ^ 

Data  for  these  studies  were  recorded  using  WARF  during  the 
months  of  March  and  November,  1974.  Observations  with  the  full 
WARF  aperture  were  taken  and  included  digitized  recordings 
from  eight  subarray  outputs  with  each  subarray  consisting  of  a 
filled  32  element  linear  array.  Both  eastward  and  westward 
looking  tests  were  included  in  the  data.  Computer  processing 
of  these  data  was  carried  out  using  the  CDC  6400  computing 
facility  at  the  University  of  Colorado  as  well  as  several  small- 
er computers  within  the  Electrical  Engineering  Department  at 
Colorado . 

The  experimental  WARF  configuration  used  to  collect  the 

r 2 t l 

data  was  identical  with  that  reported  previously1  ' . A 

sawtooth  sweep-frequency  continuous-wave  transmit  signal  similar 


to  that  shown  in  Fig.  1 was  employed  for  all  tests.  The  char- 
acteristics of  this  waveform  have  been  reported  previously ^ ' 4 ^ . 
Briefly,  at  each  receiving  subarray,  the  signal  was  deramped 
using  a local  oscillator  signal  similar  to  that  shown  in  Fig. 

1 with  a reset  time  delayed  with  respect  to  the  transmitted 
waveform.  The  delay  was  chosen  such  that  signals  from  the 
desired  range  fell  within  the  780  Hz  receiver  bandwidth  while 
those  from  other  ranges  were  rejected. 

An  analog--to-digital  converter  operating  at  a sampling 
rate  of  1.920  kHz  (32  samples  per  sweep)  was  used  to  digitize 
the  deramped  signals  after  they  had  been  converted  to  baseband 
frequencies.  Digitizing  proceeded  continuously  in  time  with  all 
samples  being  recorded  on  magnetic  tape  without  interruption. 

A multiplexed  sample  and  hold  system  provided  simultaneous 
sampling  of  the  eight  subarray  received  outputs.  Thus,  the  net 
digitizing  rate  was  8X1920=15.36  kHz.  for  the  entire  system. 
Beamformed  sampled  output  signals  were  generated  at  a 

later  time  from  these  recordings  using  the  structure  shown  in 

th  th 

Fig.  2.  If  x^(k)  is  the  k sample  from  the  i subarray,  the 
beamformed  output  signal  y(k)  is  given  by 


y(k)  = il1  W..  (k)x.  (k-j)  , 


(1) 


t h 

where  (k)  is  the  weighting  coefficient  at  the  i subarray 
which  multiplies  the  signal  after  j units  of  delay.  In  general, 
these  coefficients  are  time  varying  and  determined  by  an  adap- 
tive algorithm.  Thus,  Vh  ^ (k)  represents  the  coefficient  after 
k samples  of  data  have  been  received  from  each  subarray  and  a 
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corresponding  k adaptations  of  each  coefficient  have  taken  place. 
Adaptations  normally  are  carried  out  prior  to  the  output  compu- 
tation given  in  (1).  For  example,  in  the  case  of  the  P-vector 
algorithm,  the  adaptation  procedure  is 

wij  (k+D  = Wij  (k)+u[pi;.-y  (k)  ]Xi  (k-j)  (2) 

where  p^_.  is  the  steering  vector  coefficient.  A complete  de- 
scription of  this  algorithm  and  its  convergence  properties  are 
presented  in  Reference  3. 

As  shown  in  Figure  2 and  Eq.  (1),  a total  of  five  coeffi- 
cients per  subarray  were  used  in  the  processor.  This  number  was 
held  constant  throughout  the  results  presented  herein,  and  was 

selected  on  the  basis  of  the  recommendations  presented  in  a 
r 3 1 

previous  study  1 J . 

Performance  evaluation  of  adaptive  beamformer  performance 
was  conducted  on  the  basis  of  range-doppler  processing  of  the 
output  signal  y(k).  As  previously  discussed [ 3 1 , a two-dimensional 
digital  Fourier  transform  of  64  successive  output  sweeps  (2048 
samples  total)  results  in  a display  with  16  resolvable  range 
cells  and  64  doppler  cells.  The  resulting  display  represents 
range-doppler  information  taken  over  the  past  one  second  of  real 
time.  Figure  3 shows  a typical  range-doppler  map  derived  in 
this  manner.  The  plot  has  been  normalized  to  provide  an  average 
noise  fJoor  level  which  is  15  db  below  the  dynamic  range  of  the 
display.  Signals  which  exceed  this  range,  such  as  the  strong 
clutcer  returns  discernible  at  zero  doppler  and  all  ranges, 
are  clipped  at  the  maximum  plot  value.  The  data  for  Fig.  3 were 
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taken  on  November  8,  1974  using  eastward  transmissions  at  15.37 
MHz.  A strong  target  having  a negative  doppler  frequency  of 
about  -10Hz  can  be  seen  centered  in  range  cell  four.  A conven- 
tional beamforming  system  was  used  to  generate  this  result.  Co: 
ventional  processing  consists  of  using  fixed  weighting  coeffi- 
cients in  the  structure  shown  in  Fig.  2 with  the  following 

[4i 

Dolph  taper  coefficients: 

Wl,0  - W8,0  = °-355 

Vo  = U7,0  = °'562 
Vo  = Vo  “ °'841 
Vo  = Vo  = 1-°°° 

W.  . = 0 all  i,  j > 0. 

1,3 

Plots  of  this  type  were  computed  for  both  conventional  and  adap- 
tive beamforming  processors  under  a variety  of  signal,  noise 
and  interference  conditions  and  are  summarized  in  the  sections 
following . 


III.  COEFFICIENT  TIME-VARIATION  STUDIES 


The  adaptive  procedures  of  interest  in  this  study  are 
termed  "time-domain"  methods  to  indicate  that  the  beamforming 
coefficients  in  the  processor  are  changed  as  each  new  data 
sample  is  received.  This  is  in  contrast  with  matrix-inversion 
adaptive  methods  which  are  termed  "block-average"  techniques. 

In  these  latter  methods,  correlation  properties  of  the  received 
data  are  obtained  through  time  averages  which  are  generally 
taken  over  one  dwell — i.e.,  over  one  second  for  the  plot  shown 
in  Fig.  3.  The  resulting  averages  are  used  to  construct  a 
data  autocorrelation  matrix  Rxx  which  is  then  inverted  and  mul- 
tiplied by  the  steering  vector  Pxd  to  generate  a set  of  optimal 
coefficients  W*  (see  Eq.  (15)  in  [3]). 

w*  “ Rxxpxa  (3) 


This  set  of  coefficients  is  used  to  combine  the  subarray 
signals  observed  over  the  dwell  of  interest.  Thus,  the  output 
y(k)  in  Fig.  2 for  the  2048  samples  of  the  dwell  is  computed 
as 


y (k) 


. E.  ■ £n  w.  .x.  (k-j ) 
i=l  j=0  131 


(4) 


where  the  w^  are  determined  by  (3).  The  procedure  is  then 
repeated  at  the  next  dwell  with  a new  set  of  correlations  Rxx 
and  resulting  weights  W*.  The  procedure  is  adaptive  in  the 
sense  that  a different  set  of  weights  is  used  to  process  each 
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a 


dwell  and  the  iteration  interval  of  adaptation  is  equal  to  one 
dwell  (2048  samples) . 


This  is  in  marked  contrast  to  the  P-vector  algorithm 
method  of  adaptation.  In  matrix  form  (see  Eq.  (19)  in  [3]), 
the  algorithm  (2)  is  given  by 


where 


W (k+1)  = W(k)+u[Pxd-y  (k)  ]x(k) 


Wl,0(k) 

W2,0(k) 


w(k)  = 


W8,0(k) 
Wl,  1 {k) 


W8,4  (k) 


(5) 


(6) 


x1  (k) 
x2  (k) 


X(k) 


xg  (k) 
x^ (k-1 ) 


Xg (k-4 ) 


(7) 


(31 

It  can  be  shown  that  under  conditions  of  a stationary  environ- 
ment, the  P-vector  algorithm  (5)  and  the  matrix  inversion  method 
will  produce,  on  the  average,  identical  coefficient  vectors, 
i . e . 
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lim  E [W (k ) ] 
k 


RHS [W* ] 


(8) 


where  E [ - ] denotes  expectation. 

The  primary  advantage  of  adaptive  array  processing  methods 
in  HF  applications  is  their  capability  to  adapt  to  time-variations 
in  HF  signals  which  are  caused  by  ionospheric  fluctuations.  In 
effect,  the  processors  act  as  time-varying  linear  filters  and 
provide  interference  immunity  by  continually  changing  the  array 
pattern  so  as  to  minimize  the  effects  of  noise  and  interference 
at  the  beamformed  output.  The  effectiveness  of  an  adaptive 
beamformer  in  providing  null  tracking  is  clearly  dependent  upon 
the  time  scale  of  the  ionospheric  changes  with  respect  to  the 
adaptation  update  interval. 

Figures  4 through  8 show  the  range-doppler  maps  which  were 
computed  for  data  collected  on  March  15,  1974.  The  transmitted 
signals  for  these  data  were  directed  eastward  at  a center  fre- 
quency of  12.37  MHz.  Plots  for  both  conventional  and  P-vector 
processing  are  presented.  A proportionality  constant  of  a = 0.1 
(see  Eq.  (27),  Ref.  3)  was  used  for  these  examples.  The  target 
of  interest  is  visible  on  range  cells  6 and  7 at  a doppler 
frequency  of  about  -15  Hz.  The  strong  bands  of  energy  observed 
at  -21,  -9,  +7,  and  +23  Hz  doppler  frequencies  are  interference 
lines.  This  test  was  intentionally  conducted  under  conditions 
of  strong  HF  interference  to  illustrate  the  capabilities  of  the 
adaptive  processor  and  the  advantages  are  readily  apparent  in 
these  figures. 

Coefficient  time-variation  studies  were  carried  out  for 
the  adaptive  coefficients  used  to  generate  the  plots  shown  in 
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a)  Conventional  processing 
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DOPPLER 


b)  P-vector  adaptive  processing,  a=0.1 


Fig.  7.  Range-doppler  map  for  dwell  no.  4,  March  15,  1974 


Figs.  4b)  to  8b).  As  the  processor  adapted  on  these  data,  the 
mean  and  standard  deviation  of  each  coefficient  in  the  filter 
was  computed  at  intervals  of  128  samples  or  4 sweeps.  Since 
each  dwell  represents  an  average  over  64  successive  sweeps,  a 
total  of  16  mean  and  standard  deviation  values  were  computed 
for  each  coefficient  per  dwell.  The  results  of  this  calcula- 
tion are  presented  in  Figs.  9 and  10.  Figure  9 presents  the 
values  observed  for  the  five  coefficients  used  to  process  the 
output  of  subarray  #1  and  Fig.  10  illustrates  the  values  ob- 
served at  the  first  coefficient  or  elements  #1  through  #5. 
Standard  deviation  bars  are  shown  less  frequently  than  the  mean 
values  for  purposes  of  clarity.  Although  graphs  are  presented 
ior  only  9 of  the  40  weights  in  the  processor,  similar  curves 
were  observed  for  other  coefficients  as  well  as  for  other  exper- 
imental test  configurations. 

Computations  were  also  carried  out  to  determine  the  aver- 
age behavior  of  the  coefficient  values  over  one  complete  dwell. 
Table  I summarizes  representative  results  obtained  in  this 
manner.  In  each  case,  the  mean  and  standard  deviations  were 
carried  out  using  the  following  averaging  procedure: 


1 

16 

1 

16 

T,  i 

m=l 

Mij  (m) 

(9) 

16 

1 

16 

E, 

m=l 

°ij <ra) 

(10) 

where  (m)  and  an ^ (m)  are  the  coefficient  128  point  sample 
mean  and  standard  deviations  within  the  dwell  of  interest,  as 
shown  in  Figs.  9 and  10.  For  example,  with  m=l  and  using  the 
first  dwell, 


Table  I 


Mean  and  Standard  Deviation  Coefficient  Values 
Averaged  Over  One  Dwell 


coefficient 

dwell  no. 

mean 

std  dev 

% dev 

W10 

1 

-9.24 

.069 

.75 

2 

-8.04 

.093 

1.15 

3 

-4.25 

.093 

2.19 

4 

2.32 

.165 

7.10 

W12 

1 

6.16 

.078 

1.27 

2 

1.87 

.106 

5.63 

3 

4.34 

.066 

1.53 

4 

2.92 

.104 

3.57 

W 

14 

1 

-9.65 

.066 

.68 

2 

-7.10 

.069 

.97 

3 

-9.88 

.097 

.98 

4 

-6.85 

.089 

1.31 

W 

50 

1 

-4.88 

.074 

1.52 

2 

-6.65 

. 108 

1.62 

3 

-10.04 

.093 

.93 

4 

-7.13 

.173 

2.42 

23 
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It  should  be  noted  that  o\..  in  (10)  is  not  equal  to  the  standard 
deviation  of  the  coefficient  over  one  complete  dwell.  This  par- 
ticular method  of  averaging  is  used  here  to  eliminate  variations 
caused  by  the  time  varying  means  shown  in  Figs.  9 and  10. 

An  "average  coefficient  variation"  y similar  to  that  de- 
fined in  Equation  (35)  of  Reference  3 was  computed  for  the 
four  dwells  shown  in  Figs.  9 and  10.  The  results,  expressed 
as  a percentage,  are  given  in  Table  II. 


Table  II 

Average  Coefficient  Variation  y 


$ 


s 


Dwell  No. 
1 
2 
3 


Y 

1.06 

2.34 

1.41 


3.54 


Several  observations  are  apparent  from  these  results.  The 
large  scale  mean  weight  variations  shown  in  Figs.  9 and  10  can 
only  be  attributed  to  statistical  changes  in  the  received  data. 
Variations  caused  by  use  of  the  time  domain  P— vector  algorithm 
produce  measured  fluctuations  which  are  the  order  of  a few  per- 
cent. These  values  are  consistent  with  theoretical  values  pre- 

r 3 l 

viously  published  for  this  algorithm1  under  conditions  of 
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stationary  input  statistics  and  an  a value  of  0.1.  In  effect, 
the  time  variation  of  the  input  statistics  is  sufficiently  slow 
that  adaptive  tracking  is  excellent.  The  agreement  between 
the  y values  in  Table  II  above  and  that  published  in  Ref.  3 is 
further  evidence  of  this  fact.  Inspection  of  Figs.  9 and  10 
also  reveals  that  while  the  input  time  variation  is  slow  with 
respect  to  the  adaptive  time  constant,  it  is  of  the  same  time 
scale  as  the  dwell  interval.  Several  cases  of  marked  variations 
within  one  second  are  apparent  in  these  results.  For  the  case 
of  a two  second  dwell,  such  has  been  proposed  for  HF  backscatter 
radar  applications,  this  variation  would  be  even  more  notice- 
able. Further  discussion  and  conclusions  based  on  these  obser- 
vations are  presented  in  Section  VII  below. 


IV.  CLUTTER  SUPPRESSION  STUDIES 


In  an  HF  backscatter  radar  system,  signals  returned  from 
stationary  targets  such  as  ground  backscatter  appear  in  the 
processed  range-doppler  map  as  being  spread  over  all  range  bins 
and  centered  at  0Hz  doppler  with  a doppler  spread  of  about  1 Hz 
due  to  ionospheric  variations.  Figs.  4 through  8 show  several 
examples  of  ground-backscatter  returns,  generally  termed  clutter 
returns.  One  common  method  of  removing  these  signals  from  the 
display  is  through  the  use  of  a moving  target  indication  (MTI ) 
filter  at  the  beamformer  output.  The  general  approach  in  the 
use  of  such  filters  is  to  combine  signals  from  J adjacent  sweeps 
using  coefficients  c^.  Thus,  if  y^  (k)  is  the  kth  time  sample  of 

4-  V\ 

the  beamformed  output  during  the  j sweep,  with  index  k having 
a zero  value  at  the  start  of  each  sweep,  an  MTI  output  signal 
Zj  (k)  for  sweep  j is  formed  as 

Vk)  = /ocJYj-i(k)  (13) 

The  simplest  case  occurs  when  Cq=+  1/2  and  c^=-  1/2,  the 
two-pulse  MTI  filter,  and  adjacent  sweeps  are  subtracted.  Sta- 
tionary targets  are  eliminated  due  to  the  fact  that  they  have  the 
same  received  waveform  in  each  sweep.  More  complex  filtering 
operations  are  best  described  in  terms  of  the  effect  of  the 
operation  in  Eq.(13)  on  the  resulting  range-doppler  display. 

Since  the  c^  coefficients  are  not  functions  of  k,  they  must  have 
a constant  effect  in  the  range  domain.  Inspection  of  (13) , 
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however,  indicates  that  these  coefficients  perform  a convolution 
filtering  operation  across  the  sweeps,  for  each  specific  value 
of  k.  Since  doppler  information  is  obtained  by  a Fourier 
Transform  across  the  sweeps,  the  effect  of  this  filter  in  the 
doppler  domain  is  to  multiply  the  doppler  spectrum  by  a function 
o(wd)  given  by 


n 1 1 \ \ — JV  ^ j2TUGd£  /to  A 

c(V  ~ zt0  C£e  0 (14) 

where  wd  is  the  doppler  frequency  of  interest  and  uQ  is  deter- 
mined by  the  sweep  repetition  rate.  For  the  60  Hz  repetition 
used  in  the  present  study,  w0  = 2r  X 60  radians/sec.  In  summary, 
the  effect  of  the  MTI  operation  defined  in  (13)  is  to  multiply 
the  doppler  spectrum  for  each  range  bin  by  the  filtering  func- 
tion defined  in  Eq.  (14) . 

Thus,  for  the  two-pulse  MTI  filter,  the  doppler  filter  is 

c(wd)  = i - ie"j27Ta)d/w0  (15) 

I c (wd)  | = | sin  *wd/u>0|  (16) 

The  magnitude  function  for  this  filter  is  shown  in  Fig.  11. 

Note  that  signals  received  at  zero  doppler  are  perfectly  atten- 
uated ' c (0)  | = 0 while  those  received  at  the  doppler  fold  over 
frequencies  +a3Q/,2are  unaffected  by  the  MTI  operation  |c(iu>  /2)|  = 1. 
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Fig.  11.  Doppler  spectrum  equivalent  filter  for  two-pulse 
MTI  clutter  suppression 

The  problem  of  designing  clutter  rejection  filters  with  a 
desired  magnitude  function  c (oj^)  is  a classical  problem  in 
digital  filter  design  and  many  references  on  the  subject  are 
available. L ' J In  general,  larger  values  of  J provide  greater 
flexibility  in  this  design  but  also  impose  a longer  build-up 
time  in  the  filter.  For  example,  a 32  pulse  MTI  filter  would 
require  32  sweeps  or  about  1/2  second  to  achieve  cancellation. 

The  MTI  filters  used  in  the  present  study  were  three  pulse 
with  coefficients  given  by 

CQ  = 1/4 
C1  = - 1/2 
C2  = 1/4 

The  equivalent  doppler  function  is  then 

ctaj  = !-  se-j2™d/u>0  + je^4”^0  (17) 

d 
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I C (wd)  I = I sin  7T0J d/a) 0 1 


(18) 


and  is  plotted  in  Fig.  12. 


Fig.  12.  Doppler  spectrum  equivalent  filter  for  three- 
pulse  MTI  clutter  suppression 

Clutter  Suppression  in  Adaptive  Beamforming 

The  primary  advantage  of  MTI  clutter  suppression  in  conven- 
tional radar  systems  is  the  available  increase  in  dynamic  range. 
Removal  of  the  clutter  energy  prior  to  analog-to-digi tal  conver- 
sion allows  more  efficient  use  of  the  digitizer  in  that  the  gain 
level  is  controlled  by  desired  signals  and/or  interference,  which 
may  be  40  db  below  the  clutter.  Thus,  a larger  number  of  bits 
can  be  devoted  to  conversion  of  the  signal  and  interference 
components  than  is  the  case  when  the  gain  level  is  controlled 
by  the  clutter  energy. 


In  the  case  of  adaptive  beamforming,  there  is  an  additional 
possible  advantage  of  MTI  processing.  Figure  13  shows  a proposed 
configuration  for  clutter  suppression  in  an  adaptive  beamf ormer . 
Note  that  because  the  MTI  filters  are  placed  after  A/D  conver- 
sion, no  advantage  in  dynamic  range  can  be  achieved  A digital 
three-pulse  MTI  processor  was  placed  at  the  sampled  output  of 
each  of  the  eight  subarrays.  Thus,  the  signal  which  reached  the 
adaptive  tapped-delay  line  filter  attached  to  each  subarray  had 
little  clutter  energy  present.  The  advantage  of  this  implementa- 
tion is  that  the  adaptive  degrees  of  freedom  in  the  beamformer 
can  be  applied  directly  to  interference  suppression  and  need  not 
be  taken  up  with  the  elimination  of  clutter.  A series  of  calcu- 
lations were  carried  out  to  test  this  hypothesis.  The  data  used 
were  those  shown  in  Figs.  4-8.  Identical  data  for  each  dwell 
were  processed  with  and  without  the  MTI  filters.  In  both  cases, 
P-vector  adaptation  was  employed  with  a = 0.1.  The  range-doppler 
comparisons  for  the  five  dwells  of  interest  are  presented  in 
Figs.  14-18. 

In  order  to  better  analyze  the  behavior  of  these  MTI  filters 
signal  to  noise  ratio  and  signal  to  interference  plots  were 
derived  from  these  data.  The  measurements  were  carried  out, 
respectively,  by  plotting  received  amplitude  level  as  a function 
of  range,  with  doppler  fixed  at  the  target  doppler  and  secondly 
as  a function  of  doppler  with  range  fixed  at  the  target  range. 

Two  such  plots,  computed  for  dwell  no.  3,  are  shown  in  Figs.  19 
and  20.  In  effect,  the  plots  represent  two  perpendicular  cross- 
sections  taken  through  the  range  doppler  map  presented  in  Fig.  16 


Signal  strength  dB 
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\ ith  the  sections  crossing  on  the  target  maximum.  Finally, 
signal-to-noise  (SNR)  and  signal-to-interference  (SIR)  graphs 
were  obtained  by  extrapolating  the  noise  floor  through  the  tar- 
get signal  and  by  measuring  the  peak  interference  level,  re- 
spectively. For  the  case  of  MTI  processing,  the  interference 
levels  were  corrected  to  account  for  the  3-pulse  MTI  filtering 
function  shown  in  Fig.  12.  The  results,  plotted  as  a function 
of  dwell  number,  are  given  in  Figs.  21  and  22.  Further  discus- 
sion of  these  results  is  presented  in  Section  VII  following. 
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V.  APERTURE  TRADEOFF  STUDIES 


It  is  well  known  that  in  conventional  antenna  array  pro- 
cessing of  planar  signals  received  in  an  isotropic  noise  back- 
ground, increased  signal-to-noise  ratio  can  be  achieved  by  in- 
creasing the  receiving  array  aperture  and  maintaining  interelement 
spacing.  The  amount  of  the  increase  is  approximately  proportional 
to  N,  the  number  of  elements  in  the  total  array.  Clearly,  the 
increase  in  received  SNR  cannot  proceed  indefinitely  as  elements 
are  added  to  increase  the  receiver  aperture.  The  upper  bound  on 
output  SNR  is  limited  by  the  spatial  coherence  length  of  the 
signal  waveform.  It  has  been  demonstrated  [9]  that  at  HF , cohe- 
rence lengths  of  15  km  or  greater  can  be  realized  for  one-hop 
ionospheric  propagation  provided  that  mode  separation  methods 
are  available.  For  the  FM-CW  transmitter  format  used  in  the  pre- 
sent studies,  mode  separation  is  achieved  by  frequt  icy  filtering 
and  large  receiving  array  apertures  can  be  utilized. 

Because  the  noise  background  at  HF  is  highly  non-isotropic, 
it  is  of  interest  to  study  aperture/output  SNR  tradeoffs  for  the 
case  of  signals  received  in  the  presence  of  directional  interfer- 
ence. It  is,  in  fact,  the  non-isotropic  nature  of  the  interference 
environment  which  allows  an  adaptive  beamformer  to  provide  increased 
output  SNR  over  that  provided  by  conventional  beamforming.  Under 
conditions  of  high  levels  of  received  interference,  it  is  conceiv- 
able that  an  adaptive  beamformer  would  have  an  output  SNR  which  is 
less  sensitive  to  the  aperture  length  than  a conventional  beamformer. 
To  illustrate  this  property,  consider  the  simple  case  of  a linear, 
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equally  spaced  array  of  N isotropic  elements  spaced  by  one-half 

wavelength.  The  sample  of  the  signal  received  by  the  i^^1 

I 

- 

element  is  modeled  as  x. (k) 

\ 

x i (k)  - n (k)  + s(k)  + d (k  - iA)  (19) 

i- 


where  rw  (k)  is  white  isotropic  noise,  assumed  to  be  independent 
from  element  to  element,  s (k)  is  the  desired  signal  (sinusoidal) 
arriving  on  boresight  and  d(k  - iA)  is  a sinusoidal  interference, 
at  the  same  frequency  as  s (k)  but  incident  at  an  angle  0 which 
provides  a relative  time  delay  of  A seconds  between  adjacent  ele- 
ments. The  assumed  power  levels  of  n^ (k) , s (k)  and  d (k  - iA)  are 
°n'  °s'  and  ad'  resPectively ' at  each  of  the  N receiving  elements. 
We  consider  the  case  of  a beamformer  which  uses  amplitude 

shading  only  at  each  element.  Thus,  if  w^  is  the  weight  at  the 
t h 

i element,  the  beamformed  output  signal  y(k)  is 
n-1 

y (k ) = l w.  n . (k)  + w . s(k)  + w.  d(k-iA).  (20) 

i=0  1 


Assuming  independence  of  noise,  signal  and  interference,  the  out 
put  power  E[y2(k)]  may  be  expressed  as  a sum  of  a white  noise 
component,  desired  signal  component  and  interference  component. 


E [y 2 (k)  ] = o 
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(21) 


The  conventional  beamformer  which  is  optimal  for  the  case  of  iso- 
tropic noise  only  is  achieved  by  setting  all  weights  to  the  same 
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value.  Since  SNR  values  are  unaffected  by  scale  changes  in  the 

n--l 

weights,  a normalizing  value  providing  E w.  = 1 was  used  for 

! i=0  1 

convenience.  Thus,  with  w±  = ±,  Eq.  (21)  for  the  conventional 
processor  becomes 


, o2  o2  n-1  ... 

conv  E [y 2 (k)  ] = + o2  + -4  E e-:,Al 

N s N . . 

i = 0 


The  ratios  of  signal  to  isotropic  noise  (SNR) , signal  to  inter- 
ference (SIR)  and  signal  to  total  noise  (STR)  observed  at  the  out- 
put of  the  conventional  processor  are  then  given  by 


SNR  = 


SIR  = 


n-1 

°d  * 

i=0 
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STR  = 
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^ representative  plot  of  these  ratios,  as  a function  of  N is  shown 
in  Figure  23  for  the  case  of  o2  = 1.0,  o2  = 0.1,  o2  = 2.0,  and 
A = 30°.  Note  that  the  SNR  is  linear  in  N and  that  the  SIR  goes  to 
infinity  at  N =12.  This  corresponds  to  an  aperture  such  that  a 
null  of  the  conventional  pattern  lies  in  the  direction  of  the  inter- 
ference component. 

An  adaptive  beamformer  for  the  same  problem  was  also  imple- 
mented. In  this  case,  the  amplitude  weights  w^  used  for  each  element 
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output  were  chosen  using  a minimum  mean-square  error  criterion. 
Thus,  the  vector  of  weight  coefficients  W was  computed  using  the 
optimal  solution  given  by  Eq.  (3)  in  Section  III  above.  For  the 
simple  model  presently  under  consideration,  the  autocorrelation 
matrix  and  steering  vector  P are  given  by 

RYY  = o2I  + o211T  + a2 A (24) 

xx  n s d 

P = a2l  (25) 

s 

where  1 is  a column  vector  of  one's  and  A is  an  NXN  matrix  in 
which  the  (i,j)  element  is  cos[A(i-j)].  The  optimal  value  of  W 
was  computed  for  N ranging  from  1 to  20.  The  appropriate  ratios 
SNR,  SIR,  and  STR  were  then  computed  using  the  individual  terms 
computed  in  Eq.  (21).  The  results,  plotted  as  a function  of  the 
number  of  elements  N,  are  presented  in  Fig.  24. 

A comparison  of  Figs.  23  and  24  reveals  several  pertinent 
properties  relating  to  the  importance  of  aperture  in  conventional 
and  adaptive  beamformers.  Although  these  results  were  generated 
using  amplitude  taper  only,  they  are  representative  of  the  effects 
which  are  observed  using  either  tapped-delay-line  processing  or 
gain-phase  processing.  The  first  observation  is  that  the  ratio  of 
desired  signal  power  to  undesired  signal  power,  STR,  observed  at 
the  output  of  a conventional  beamformer  does  not  increase  linearly 
with  increasing  aperture  when  the  noise  is  isotropic.  While  the 
trend  is  generally  such  that  increased  aperture  provides  increased 
STR,  there  are  cases  where  a reduction  in  aperture  increases  this 
ratio  - i.e.,  for  N in  the  range  12  to  17.  The  second  observation 


44 


is  that  adaptive  processing  provides  a monotonic  increase  in  STR 
with  increasing  aperture  and  that  the  value  of  STR  achieved  very 
nearly  saturates  after  a certain  minimum  aperture  has  been  achieved. 
The  reason  for  the  monotonic  behavior  is  simply  that  as  more  ele- 
ments are  added,  additional  degrees  of  freedom  become  available 
for  cancelling  the  directional  interference.  The  saturation  effect 
is  caused  by  the  fact  that  once  a certain  minimum  number  of  degrees 
of  freedom  have  been  obtained,  further  reduction  in  interference 
power  can  be  achieved  only  at  the  expense  of  increased  isotropic 
power  at  the  outout  - i.e.,  compare  the  first  and  third  terms  in 
Eq.  (21).  This  saturation  effect  has  been  observed  elsewhere  [8] 
and  has  been  shown  to  be  a function  of  the  absolute  isotropic  noise 
power  level.  That  is,  when  simulations  are  carried  out  using  values 
of  lower  than  that  selected  for  the  present  example,  the  satura- 
tion effect  occurs  at  a higher  value  of  N and  the  resulting  STR 
value  is  also  greater.  One  factor  clearly  demonstrated  by  Fig.  24 
is  that  under  some  conditions,  a reduction  in  aperture  may  have 
virtually  no  effect  on  output  STR  for  an  adaptive  processor. 

Experimental  Aperture  Tradeoff  Studies 

In  order  to  investigate  the  effect  of  aperture  changes  on  the 
performance  of  an  adaptive  processor  operating  in  an  actual  HF  sig- 
nal environment,  the  data  collected  on  March  15,  1974  and  previously 
shown  in  Figs.  4 through  8,  was  processed  at  less  than  full  aper- 
ture. In  each  case,  P-vector  adaptation  was  used  with  a = 0.1  and 
5 adaptive  coefficients  per  subarray  output.  Figure  25  shows 
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results  obtained  during  dwell  No.  4 with  3-pulse  MTI  in  use.  The 
three  plots  represent  processing  of  the  same  data  values  with  two 
adjacent  subarrays  (Fig.  25a),  four  adjacent  subarrays  (Fig.  25b), 
and  the  full  array  of  eight  subarrays  (Fig.  25c). 

Comparison  of  Figs.  25b  and  c indicates  that  reduction  of  the 
aperture  size  by  a factor  of  two  had  little  detrimental  effect  on 
array  performances.  Further  reduction  by  an  additional  factor  of 
two,  however,  produced  significantly  inferior  results.  The  effects 
of  aperture  were  further  studied  by  varying  the  aperture  size  while 
maintaining  the  number  of  subarrays  at  two.  Figure  26  shows  re- 
sults obtained  using  two  adjacent  subarrays,  two  at  half  aperture 
and  two  at  the  full  aperture.  The  two  at  half  aperture  illustrate 
a clear  advantage  over  any  other  aperture.  This  effect  was  not  re- 
stricted to  adaptive  processing,  as  illustrated  in  Fig.  27.  In 
this  case,  the  same  data  were  processed  using  a conventional  beam- 
former.  Again,  the  superiority  of  the  half  aperture  configuration 
is  apparent.  The  effect  was  also  not  observed  on  only  one  dwell 
but  occurred  frequently  throughout  the  data  record.  Figure  28 
illustrates  a similar  effect  obtained  using  the  adaptive  processor 
on  dwell  no.  3. 

One  possible  explanation  for  this  effect  is  the  presence  of 
a fading  null  across  the  array.  Since  the  scale  size  of  HF  fading 
structures  is  known  to  be  smaller  than  the  Los  Banos  aperture,  it 
is  possible  that  several  of  the  subarrays  were  resident  in  a fading 
null  during  these  tests.  It  is  clear  that  the  problem  of  aperture 
tradeoffs  in  the  HF  environment  is  not  a simple  one  and  further 
study  of  these  effects  is  warranted.  The  phenomenon  illustrated  in 


47 


Figs.  25-28  are  further  evidence  of  the  advisability  of  using 
adaptive  processing  at  HF.  Since  such  effects  are  most  certainly 
time  varying,  one  advantage  offered  by  adaptation  is  that  of  placing 
greater  weight  on  the  high  SNR  subarray  outputs. 


VI.  MINIMUM  MEAN-SQUARE  ERROR  ADAPTATION  WITH  LINEAR  CONSTRAINTS 


The  adaptive  beamforming  method  used  to  process  data  in  this 
and  previous  ' reports  has  been  referred  to  as  the  P-vector 
adaptation  algorithm.  A complete  description  of  this  algorithm, 
including  its  derivation  and  convergence  properties,  has  been 
presented  elsewhere.  Briefly,  the  procedure  is  designed  to 

minimize  the  mean-square  error  between  a desired  signal  d(k), 
which  is  presumed  to  be  incident  on  the  array  from  a known  direc- 
tion. For  the  FM/CW  backscatter  radar  case,  the  desired  response 
is  assumed  to  be  sinusoidal  with  a frequency  determined  by 
the  range  of  interest.  Thus,  the  desired  response  is 

d (k)  = cos  (wQk)  (26) 

and  the  signal  observed  at  the  sampled  output  of  the  ith  subarray 
is  assumed  to  be  of  the  form 

Xi(k)  = Aid{k-ii)  + ni  (k)  (27) 

where  n^  (k)  is  the  sum  of  all  noise  and  interference  components 
t h 

at  the  i subarray  output.  Under  these  assumptions,  the  beam- 
formed  output  signal,  y(k)  [Eq.(l)J,  can  be  expressed  as  a sum  of 
desired  and  noise  components, 

Y(k)  = Yd(k)  + Y (k ) (28) 

The  P-vector  algorithm,  Eq ' s (2)  and  (5),  is  then  designed  to 
adapt  toward  a set  of  weighting  coefficients  which  minimize  the 
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mean-square  difference  E(e2(k)]  between  d(k)  and  Y(k),  i.e. 

E [e  2 (k)  ] - E [{ d (k)  - Y (k ) } 2 ] (29) 

The  optimum  vector  of  coefficients  W*  which  minimizes  (29) 
is  given  by  Eq.  (3)  and  is  referred  to  as  the  constraint-free, 
minimum  mean-square  error  beamforming  solution.  One  significant 
property  of  this  solution  is  that  the  resulting  array  processor 
may  provide  frequency-domain  filtering  of  signals  received  from 
the  desired  direction.  An  example  of  this  filtering  is  illus- 
trated in  Figure  29.  This  two-dimensional  plot  shows  the  array 
pattern  response  at  16  resolvable  azimuths,  labeled  0 through  15, 
as  a function  oi  frequency  within  the  receiver  bandwidth,  0 to 
960  Hz.  The  plot  was  obtained  using  a two-dimensional  Fourier 
Transform  of  coefficients  generated  by  use  of  the  P-vector 
adaptation  algorithm  on  data  recorded  September  3,  1974.  East- 
ward transmissions  were  employed  and  the  beamformer  consisted  of 
five  coefficients  per  subarray  with  a = 0.1.  MTI  processing  was 
not  in  use  for  this  test.  The  desired  direction  of  arrival  cor- 
responded to  azimuth  cell  number  8 and  the  desired  frequency  was 
480  Hz,  as  shown  by  the  arrow. 

Note  that  any  signal  incident  on  the  array  from  the  desired 
direction  and  at  the  desired  frequency  is  in  the  main  lobe  of 
the  pattern.  In  addition,  a variety  of  nulls,  extending  in  both 
frequency  and  azimuth,  are  evident.  The  location  and  nature  of 
these  nulls  is  a reflection  of  the  noise  and  interference  en- 
vironment under  which  data  were  recorded.  A variety  of  similar 

plots,  taken  under  many  other  data  configurations,  are  presented 
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in  a recent  SRI  report. 
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Inspection  of  Fig.  29  reveals  that 
on  the  cirray  from  the  desired  direction,  but  at  a frequency  other 
than  the  desired  480  Hz  value,  will  be  attenuated  with  respect  to 
the  response  at  480  Hz.  In  the  range-doppler  display,  frequency 
tapering  of  this  type  appears  as  range-dependent  filtering. 
Examination  of  the  interference  lines  in  the  adaptive  plots 
presented  earlier.  Figs.  4 through  8,  shows  further  evidence  of 
this  effect.  It  is  to  be  noted,  however,  that  the  frequency  (or 
range)  filtering  is  never  particularly  sharp  for  the  filter  para- 
meters used  in  this  study  and  will  not,  for  example,  modify  an 
interference  line  to  the  point  of  making  it  appear  target-like. 

The  sharpness  of  this  filtering  is  increased  as  the  number  of 
adaptive  coefficients  per  subarray  output  is  increased.  With 
five  taps,  the  effective  minimum  bandwidth  of  the  filter  is 
approximately  6.5  range  lines.  Extending  the  processor  to  16 
taps  per  subarray  would  reduce  the  bandwidth  to  two  range  lines 
and  might  very  well  cause  interference  to  appear  target-like  in 
the  range-doppler  display. 

It  is  possible  to  derive  optimal  beamforming  procedures 
which  are  constrained  to  have  a flat  frequency  response  in  the 
desired  look  direction  of  the  array,  regardless  of  the  number  of 
taps  used  to  process  each  subarray  output.  Such  procedures  fall 
into  a class  termed  minimum  mean-square  error  with  constraints 
and  have  been  studied  extensively  by  previous  authors.  f 11-13 J 
0f  Particular  interest  is  a recent  procedure  suggested  by  Frost ^ 5 ^ 
which  allows  time-domain  adaptation  under  linear  constraints. 

In  the  sections  following  a brief  description  of  the  optimal 
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processor  and  Frost's  algorithm  are  presented.  Comparative 
experimental  results  obtained  using  both  P-vector  and  Frost 
adaptation  on  the  same  data  are  also  given. 


Optimum  Constrained  Least-Mean  Squares  Filter  Coefficients 

The  constrained  minimum-mean  square  error  beamformer  is 
designed  to  minimize  the  total  array  output  power  subject  to 
a linear  constraint  on  the  filter  coefficients.  (Clearly,  with 
no  such  constraint,  the  optimum  solution  is  obtained  by  setting 
all  coefficients  to  zero.)  With  reference  to  the  weight  vector 
defined  in  Eq.  (6),  the  constraint  may  be  expressed  as 

CTW  = F (30) 

where  W is  the  K-L  dimensional  vector  of  weights  for  a processor 
containing  K subarrays  and  L taps  per  subarray;  C and  F are, 
respectively,  a J by  KL  matrix  and  a J-dimensional  vector  of 
known  constants  which  determine  the  constraints.  In  effect,  Eq . 
(30)  is  a specification  of  J < KL  simultaneous  constraint  equa- 
tions on  the  weight  vector. 

For  array  processing  with  a frequency  independent  main  lobe, 
the  appropriate  constraint  is  that  the  sum  of  K coefficients  taken 
across  the  array  after  any  number  of  identical  delays  should  be 
zero  except  for  one  pre-speci f ied  delay  value,  which  has  a sum 
of  unity.  Mathematically,  for  the  8 subarray,  5 tap  per  subarray 
case,  this  is  expressed  as 


0 


j = 0,1, 3, 4 


w . . 
ID 


(31) 


With  this  constraint,  any  signal  which  arrives  on  boresight  is 
received  by  all  eight  subarrays  simultaneously  and  therefore 
appears  undistorted  at  the  output  when  it  reaches  tap  no.  2. 

Note  that  the  processor  requires  time  delay  beamsteering  prior 
to  processing  to  provide  an  effective  boresight  desired  look 
direction. 

The  equivalent  C and  f matrices  in  (30)  for  this  constraint 
are  given  by 


1 0 

1 0 

• • • 
1 0 
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0 0 

0 0 


(32a 
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F = 


1 
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(32b) 
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The  optimum  weight  vector  W*  which  minimizes  E[y  (k) ] subject 
to  the  constraint  (30)  is 

W*  = R^C[CTr“^C]_1F  (33) 

[12  1 

It  can  be  shown  that  this  solition  is  closely  related  to  the 

unconstrained  minimum  mean-squart  error  solution  given  by  Eq . (3) 

when  C and  F have  the  form  in  132) . In  particular,  the  uncon- 
strained solution  has  been  shown  to  have,  in  general,  a higher 
output  signal-to-noise  ratio  than  does  the  constrained  method. 

In  the  limit  of  arbitrarily  high  input  signal  to  total  noise 
ratio,  however,  the  two  processors  become  virtually  identical. 
Qualitatively,  the  lower  STR  in  the  constrained  case  is  due  to 
the  fact  that  some  of  the  adaptive  degrees  of  freedom  are  taken 
up  by  the  constraint.  The  extreme  example  of  this  would  occur 
when  the  number  of  constraints  was  set  equal  to  the  number  of 
coefficients.  Thus  in  Eq.  (30)  with  J=KL,  the  matrix  C is  square 
and,  if  non-singular,  completely  determines  the  weight  values, 
irrespective  of  the  data  received  by  the  array. 


Constrained  Least-Mean  Square  Adaptation 


r 5 1 

Frost  has  presented  a time-domain  adaptation  procedure 
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which  may  be  used  to  adapt  toward  the  optimal  solution  defined  by 
Eq.  (33).  The  method  is  similar  to  the  P-vector  algorithm  in  that 
the  processor  coefficients  are  updated  as  each  new  sample  of  data 
is  received  and  measurement  of  signal  correlation  matrices  are 
not  required.  Frost's  algorithm  has  the  additional  advantage 
that  it  is  self-correcting  and  avoids  error  accumulation  while 
maintaining  the  specified  constraint.  Roundoff  and  truncation 
errors  thus  have  a minimal  effect  upon  the  performance  of  the 
algorithm. 

For  the  general  constraint  conditions  defined  by  Eq . (30), 
Frost's  procedure  may  be  expressed  as 

W (k+1)  = Q[W(k)  -py(k)X(k)J  + F (34) 

where  Q is  a square  KL  by  KL  matrix  of  known  coefficients, 

T - 1 T 

Q = I - C (CC1 ) V (35) 

and  W(k)  and  X(k)  are  the  coefficient  and  signal  vectors  defined 
by  Eq . ' s (6)  and  (7)  above.  When  the  C matrix  has  the  particularly 
simple  form  shown  in  Eq.  (32a),  the  premultiplication  by  Q in 
(34)  can  be  expressed  in  terms  of  additions  as  follows: 


wi ,o(k+l ^ = wi  0(k)-My(k)x, 

(k)-  i[JJ1wjf0(k)-py(k)Xj(k)]+ 

(36a) 

w8  o(k+,)  = w8  0(k)-ljy(k)X8 

1 8 f , 
(k)-  q"[ j ^ , w j )0(k)_k|y (k)xj (k)  ] + -g1 

(36b) 

w,  J (k+l ) = Wj  1 (k)  — jjy  (k ) X 1 

« 8 

(k-0-  j[jx1wj  J(k)-My(k)Xj(k-l)]+  ~ 

(36c) 

w8>i4(k+D  = w8>4(k)-uy(k)X8 

1 8 fr 

(k-v-  B'[j>J1wj>i,(k)-uy(k)xj  (k-4)]+  -f 

(36d) 
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The  total  number  of  digital  operations  required  to  implement 
these  equations  on  a digital  computer  for  a 40  weight  processor 
is  95  real  multiplies  and  240  real  additions.  In  the  case  of 
P-vector  adaptation,  the  numbers  are  81  multiplies  and  120  addi- 
tions. It  should  be  noted  that  for  both  processors  these  values 
include  the  40  multiplies  and  adds  necessary  to  compute  the 
beamformed  output  y(k). 

Experimental  Comparisons 

Figures  30  through  34  show  a series  of  successive  range- 
doppler  plots  obtained  using  a one-second  integration  period. 

In  each  case,  the  same  data  were  processed  using  both  P-vector 
and  Frost  adaptation  methods  with  five  adaptive  coefficients 
per  subarray  and  an  a value  of  0.1.  The  data  were  recorded  using 
eastward  looking  transmissions  on  March  15,  1974  and  are  iden- 
tical to  that  shown  in  Figs.  4-8.  MTI  filters  were  not  in  use 
during  processing. 

A comparison  of  these  results  illustrates  that  the  received 
target  signal-to-noi se  and  signal-to-interf erence  ratio  is 
generally  greater  for  P-vector  adaptation  than  for  adaptation 
using  the  Frost  procedure.  Cross-sectional  received  signal  plots 
taken  at  the  target  doppler  frequency  as  a function  of  range  are 
shown  in  Figs.  35  and  36  for  dwell  no.  3.  Also  shown  for  pur- 
poses of  comparison  are  the  results  obtained  with  conventional 
processing . 
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a)  Frost  algorithm 


b)  P-vector  algorithm 


Fig.  33.  Constrained  minimum  mean-square  error  comparison 
for  dwell  no.  4,  March  15,  1974 
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b)  P-vector  algorithm 


Constrained  minimum  mean-square 
for  dwell  no.  5,  March  15,  1974 


amplitude  response  observed  in  velocity  cell  -17  during 


The  frequency-independent  nature  of  the  boresight  response 
of  the  Frost  procedure  is  illustrated  by  the  f requency/azimuth 
plot  presented  in  Fig.  37.  Comparison  with  the  equivalent 
P -vector  response  given  in  Fig.  29  shows  that  the  primary  differ- 
ence lies  in  the  response  observed  at  boresight  (azimuth  cell 
number  8 ) . 
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VII.  DISCUSSION  AND  CONCLUSIONS 


This  report  has  summarized  the  efforts  carried  out 
at  the  University  of  Colorado  during  the  period  March,  1974 
to  March,  1975  in  support  of  the  ARPA  HF  Adaptive  Array 
Program.  The  emphasis  at  Colorado  has  been  on  the  detailed 
analysis  of  adaptive  beamforming  characteristics,  including 
the  development  and  analysis  of  various  algorithms,  using 
a relatively  small  sample  of  field-recorded  data.  Detailed 
studies  which  test  these  procedures  over  longer  time  periods 
and  under  a variety  of  ionospheric  and  interference  condi- 
tions need  to  be  conducted  for  purposes  of  gathering  a 

statistical  history  of  adaptive  beamforming  performance. 

[ 4 1 

F<  rtunately . such  studies  have  been  conducted  at  SRI 
during  the  past  year  using  a real  time  implementation  of 
the  P-vector  algorithm  and  will  continue  over  the  next  year. 

The  results  obtained  at  Colorado  during  the  past  year 

[2  31 

have  substantiated  earlier  reported  results.  ' Speci- 
fically, it  has  been  demonstrated  that  substantial  signal- 
to-noise  and  signal-to-inter f erence  improvements  can  be 
achieved  in  HF  backscatter  radar  systems  which  employ  time- 
domain  adaptive  beamforming  methods.  It  has  been  shown 
that  interference  levels  in  the  processed  range/doppler 
maps  can  be  reduced  as  much  as  20  dB.  In  addition,  the 
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overall  noise  floor,  exclusive  of  those  areas  containing 
interference,  is  substantially  lowered.  Under  every  con- 
dition investigated  to  date,  the  adaptive  beamformer  has 
provided  performance  levels  superior  or  equal  to  that  offered 
by  conventional  processing.  No  conditions  have  been  found 
in  which  the  use  of  adaptive  processing  induced  a degrada- 
tion in  performance. 

Several  new  specific  results  were  obtained  during  the 
course  of  the  present  study  and  have  been  grouped  into 
four  areas  of  investigation: 

1 • Coefficient  Time-Variation  Studies 

It  has  been  demonstrated  that  the  time  scale  of  the 
variation  in  optimal  weights  in  an  HF  environment  is  the 
order  of  one  Second  or  less.  In  several  cases,  significant 
variability  was  noted  within  a time  span  as  short  as  one 
half  second.  It  was  further  demonstrated  that  the  P-vector 
time-domain  adaptation  procedure  can  successfully  track 
these  variations  when  used  with  a step  size  value  of  a = 0.1. 
The  standard  deviation  of  the  adaptive  weights  about  their 
time  varying  mean  values  was  found  to  be  the  order  of  2-3%. 
This  is  in  excellent  agreement  with  the  value  predicted  by 
the  theory  of  convergence  of  the  adaptive  algorithm  for 
the  case  of  stationary  input  statistics.  Thus,  although 
the  signals  observed  at  HF  are  non-s tationary , the  time 
scale  of  their  variation  is  sufficiently  slow  so  as  to 
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be  successfully  followed  by  the  P-vector  algorithm.  It 
should  also  be  noted/  however,  that  this  conclusion  holds 
only  for  a processor  consisting  of  40  adaptive  weights. 

The  theory  of  adaptive  processing  has  shown  that  adaptive 
convergence  times  are  approximately  linearly  related  to 
the  number  of  adaptive  degrees  of  freedom  in  the  system. 
Assuming  that  the  time  constant  for  the  present  system  is 
about  0.1  second,  this  places  an  upper  limit  of  about  200 
coefficients  for  an  adaptive  signal  processing  system 
operating  in  an  HF  environment.  At  this  level,  the  adaptive 
time  constant  would  be  about  0.5  seconds,  a value  insuffi- 
cient to  track  the  observed  variations. 

2 . Clutter  Suppression  Studies 

The  use  of  individual  clutter  suppression  filters  at 
the  individual  subarray  outputs  has  been  investigated.  For 
the  rather  small  data  sample  used  during  this  study  (five 
successive  one  second  dwells)  it  was  found  that  MTI  filters 
provided  neither  significant  advantage  nor  significant  dis- 
advantage in  processing.  Signal-to-noise  and  signal-to- 
interference  ratios  measured  in  the  range-doppler  maps 
varied  less  than  a few  dB  between  using  and  not  using  three 
pulse  clutter  suppression  filters.  Cases  were  found  where 
the  use  of  MTI  filters  produced  inferior  performance  and  in 
other  examples,  superior  performance  over  that  obtained  by 
adapting  in  the  presence  of  clutter.  Similar  studies  carried 
out  by  SRI ^ have  indicated  an  example  in  which  MTI 
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processing  provided  a seven  dB  advantage  (see  Fig.  12, 

Ref.  4).  This  advantage  was  observed  only  for  a two  second 
interval  out  of  a nine  second  record,  however,  and  for  the 
remainder  of  the  time,  no  significant  advantage  or  disad- 
vantage was  evident. 

It  is  difficult  to  draw  definitive  conclusions  regard- 
ing the  use  of  MTI  filters  in  an  adaptive  beamforming 
system  based  upon  the  present  evidence.  Certainly,  the 
qualitative  argument  which  states  that  MTI  filters  must 
help  due  to  the  fact  that  the  removal  of  clutter  releases 
adaptive  degrees  of  freedom  is  convincing.  At  this  point, 
it  appears  that  the  use  of  clutter  suppression  devices  is 
at  yet  unresolved  and  further  investigative  efforts  in 
this  area  are  warranted. 

3 • Aperture  Tradeoff  Studies 

The  effect  of  reduced  aperture  size  on  beamformer 
performance  has  been  studied  for  both  conventional  and 
adaptive  processing.  A simple  simulation  model  was  used 
co  demonstrate  that  in  a non-isotropic  noise  environment, 
conventional  processing  does  not  provide  increased  output 
signal-to-noise  plus  interference  ratio  as  elements  are 
added  to  the  array  to  expand  the  aperture.  This  model 
also  demonstrates  that  adaptive  beamforming  in  the  same 
environment  did  provide  a monotonically  increasing  per- 
formance with  increased  aperture  size,  but  that  a threshold 
aperture  size  existed  above  which  the  performance  increase 
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was  minimal. 


Observations  using  experimental  data  confirmed  that  a 
reduction  of  up  to  one  half  the  original  aperture  did  not 
significantly  reduce  system  performance  when  adaptive  beam- 
forming was  in  use.  It  was  further  found  that  the  rela- 
tionship between  output  SNR  and  aperture  size  for  the  case 
of  two  subarrays  with  variable  spacing  was  not  simply  related 
to  the  aperture  spanned  by  the  subarrays.  Specifically, 
processing  subarray  no's  1 and  4,  which  span  one  half  of 
the  available  Los  Banos  aperture,  produced  larger  target 
signals  in  the  range-doppler  display  than  did  the  process- 
ing of  subarrays  1 and  8 (full  aperture)  or  numbers  1 and 
2 (one-eighth  aperture) . This  result  was  observed  for  both 
conventional  and  adaptive  beamforming  methods.  One  possible 
explanation  for  these  results  is  the  prese-ce  of  a wave- 
like fading  structure  across  the  receiving  array.  Under 
such  conditions,  the  question  of  optimal  subarray  spacing 
and  location  could  be  expected  to  change  as  the  fading 
nulls  move  across  the  array.  Further  detailed  investiga- 
tions of  this  phenomena  is  warranted. 

^ • Minimum  Mean-Square  Error  Adaptation  with  Constraints 
A generalized  time-domain  adaptive  procedure  for  mini- 
mizing the  beamform^d  output  power  subject  to  a linear 
constraint  on  the  processor  coefficients  has  been  imple- 
mented and  tested  using  experimental  data.  Results  have 
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shown  that  satisfactory  convergence  rates  and  performance 
can  be  achieved  at  a computational  cost  approximately  50% 
greater  than  that  required  for  P-vector  adaptation.  The 
particular  constraint  studied  in  these  tests  provided  a 
frequency  independent  response  in  the  main  lobe  direction. 

It  was  shown  that  lower  output  signal-to-noise  and  signal- 
to-interference  ratios  are  achieved  with  this  method  than 
with  P-vector  adaptation  and  that  the  difference  may  be 
as  great  as  five  dB  in  some  range-doppler  cells. 

Suggestions  for  future  work  with  constrained  adaptive 
methods  include: 

a)  Additional  data  processing  to  provide  comparison 
with  the  P-vector  method  under  differing  ionospher- 
ic ar.d  interference  conditions. 

b)  An  investigation  into  the  utility  of  linear  con- 
straints other  than  that  of  a frequency  independent 
main  lobe. 

c)  A study  of  the  performance  of  the  constrained  pro- 
cedures under  conditions  in  which  one  or  more 
hardware  elements  fail  in  the  adaptive  processor.  It 
is  conceivable  that  procedures  designed  to  minimize 
output  power  may  very  well  place  a large  weighting 
on  a quiet  but  useless  subarray  output  produced, 

for  example,  by  failure  of  a preamplifier. 
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5 . Suggestions  for  Future  Studies 

The  following  suggestions  for  future  research  in  the 
area  of  adaptive  beamforming  at  HF  became  evident  during 
the  course  of  the  present  investigations: 

a)  A study  of  additional  time-domain  adaptive  pro- 
cedures which  are  designed  to  optimize  criteria 
other  than  minimum  mean-square  error  should  be 
conducted.  This  effort  would  include  maximizing 
output  SNR  and  the  various  methods  of  adaptive 
sidelobe  cancelling  which  have  been  suggested. 
Comparisons  with  the  existing  least-squares  tech- 
niques should  be  conducted  using  simultaneous 
processing  of  experimental  data. 

b)  A detailed  comparison  between  the  time-domain, 
single  receiver  channel,  tapped-delay-line  method 
of  processing  discussed  in  this  report  and  the 
alternative  procedure  of  quadrature  receiver 
processing  should  be  carried  out.  Optimization 
in  the  latter  case  is  generally  achieved  using 
correlation  matrix  estimation  and  inversion  methods. 
Time-domain  procedures  for  thi-  class  of  receiver 
are  also  available,  however,  and  should  be  included. 

c)  In  most  HF  installations  proposed  to  date,  the 
number  of  receiving  elements  greatly  exceeds  the 
number  of  adaptive  degrees  of  freedom  which  can  be 
utilized.  The  problem  of  optimally  distributing 
these  degrees  of  freedom  in  the  system — i.e.  the 
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subarray  geometry  problem — is  as  yet  unsolved  and 
warrants  further  study.  One  reasonable  approach 
to  this  rather  ill-defined  problem  would  be  to 
consider  only  minimum  redundancy  and  maximum 
redundancy  configurations.  The  former  case  has 
been  widely  investigated  by  radio  astronomers  but 
has  the  distinct  disadvantage  that  overall  perfor- 

I 

mance  is  greatly  deteriorated  by  the  loss  of  a 
single  element,  such  as  during  a deep  fading  null. 
Maximal  redundancy  avoids  this  problem  by  distri- 
buting the  degrees  of  freedom  uniformly  among  the 
available  elements. 

Since  any  operational  HF  radar  system  will  operate 
using  multiple  simultaneous  received  beams,  further 
research  must  be  conducted  to  determine  the  best 
hardware  configuration  for  a multibeam  adaptive 
processor.  This  study  should  include  the  possi- 
bility of  performing  separate  FFT  operations  at 
each  subarray  output,  prior  to  adaptation. 
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